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ABSTRACT 



Aims. We study the evolution of the galaxy population up to z ~ 3 as a function of its colour properties. In particular, luminosity functions and 
luminosity densities have been derived as a function of redshift for the blue/late and red/early populations. 

Methods. We use data from the GOODS-MUSIC catalogue which have typical magnitude limits zaso ^ 26 and Ks < 23.5 for most of the 
sample. About 8% of the galaxies have spectroscopic redshifts; the remaining have well calibrated photometric redshifts derived from the 
extremely wide multi-wavelength coverage in 14 bands (from the U band to the Spitzer 8/jm band). We have derived a catalogue of galaxies 
complete in rest-frame B-band, which has been divided in two subsamples according to their rest-frame U-V colour (or derived specific star 
formation rate, SSFR) properties. 

Results. We confirm a bimodality in the U-V colour and SSFR of the galaxy sample up to z ~ 3. This bimodality is used to compute the LFs of 
the blue/late and red/early subsamples. The LFs of the blue/late and total samples are well represented by steep Schechter functions evolving in 
luminosity with increasing redshifts. The volume density of the LFs of the red/early populations decreases with increasing redshift. The shape 
of the red/early LFs shows an excess of faint red dwarfs with respect to the extrapolation of a flat Schechter function and can be represented 
by the sum of two Schechter functions. Our model for galaxy formation in the hierarchical clustering scenario, which also includes external 
feedback due to a diffuse UV background, shows a general broad agreement with the LFs of both populations, the larger discrepancies being 
present at the faint end for the red population. Hints on the nature of the red dwarf population are given on the basis of their stellar mass and 
spatial distributions. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: luminosity functions 



1. INTRODUCTION 

The evolution of galaxy Luminosity Function (LF) is one of 
the main tools to study the structure evolution through the cos- 
mic time. The advent of large surveys has allowed the analysis 
of sub-samples of galaxies selected as a function of their mor- 
phological, spect roscopic or colour properties ( Strateva et al. 

2002; Wol fet all 
B~ldrv et al. 2004'; 
Blantonetal. 2005; 



2001 



2003 



Norberg et al. 



Wilkner et al. 



Hogg etal. 12004 



Ilbert et al 



Driver et al 



2002; Madgwick 



dm06; Bell e t alj|2()0 
3^iner et a lj 120051; 



ketal.ll2( 
2004'Bal 



2006; Marchesini et al. 12007 



iBaldrv et all 



2006t iCuc ciati et al."2006'; 'C irasuolo et all 



Amouts et al. 



.I l2007t 



2ooe 
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Of special interest are the studies concerning the statis- 
tical properties of galaxies selected on the basis of their in- 
trinsic colour distributio n. This d istribution appears bimodal 
up to z ~ 2-3 (Baldry et all i2004; Blantonet_al. 20051: 



iGiallongo et alJ2005l) and separates the galaxies in two popula 
tions, red early types vs. blue late types. It has been shown that 
this spectral classification is roughly consistent with the cor- 
respondent morphological classification (bulge vs. disc domi 



nated) at least at low and intermediate redshifts (IStrateva et al 
2001; Weineretal. 2005). 

This bimodal colour distribution can find a natt iral explana- 



Giallongo et al. 2005, elsewhere G05). In tion in hier archical models for galaxy formation dMenci et al 



fact these kind of studies allow us to probe the evolution of 
galaxies having different star formation histories. 



20051 12006 ) where two distinct populations arise in the colour 
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distribution based on two different star formation histories af- 
fected by the feedback effec ts produced by the SN and AGN 
activities (^ Menci etal 2006h . However the effect of environ- 
mental density on the paths of galaxy evolution can have a fun- 
damental role. In this context it is not clear whether the evo- 
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lutionary history of galaxies is originated by a nurture senario 
(galaxy properties are affected by environment through phys- 
ical mechanisms acting on galaxies) or by a nature scenario 
(the evolution is driven by the initial conditi on estab lished dur- 
ing the formation epoch of galaxies, e.g. IMateus et al.l 120071: 
Cooper et alJl2007h . 

Recent studies have estimated the shape and evolution 
of the LF of galaxies selected according to their bimodal 
colour distribution using both the large local Sloan sur- 
2004; Blanton et al. 2 0051 

and high redshi f ts (IBell et alj 



vey jBaldrv et al. 
veys at i ntermediate 



[Giallon20 et a]^hoM 
Ilbert et al. 2006). Their results show the red LF evolving 
mildly in density up to z ~ 1 with a quite flat shape at the 
faint-end, although the evaluation of the faint-end slope of the 
red LF remains an open issue especially at intermediate and 
high redshifts where the present surveys do not constr ain the 
faint slope very weU dFaber et al.ll2005HBeU et"al]|2003h . 

In G05 we studied the red and blue LFs, using the proper- 
ties of bimodality in colour and in specific star formation rate 
(SSFR), with a complete but relatively small sample of galax- 
ies selected in the rest-frame B-band from low to high redshifts. 
We showed that the bimodality extends at least up to z ~ 2.5. 
We also found that the red/early galaxies decrease in their lumi- 
nosity density by a factor ~ 5 - 6 from z ~ 0.5 to z ~ 2.5 - 3 in 
broad agreement with the hierarchical cold dark matter model. 
These results provided a first picture of the evolution of the red 
and blue LFs up to high redshifts relaying on a relatively deep 
but small sample. For a more reliable picture a wider sample 
at high redshift is clearly needed. For this purpose larger areas 
with deep near-IR imaging are required. 

Thanks to the wide area (~ 140 arcmin^) and to the deep 
near-IR observations, the GOODS-South survey provides a 
good starting point for the study of the galaxy properties at 
high redshift. In particular, the inclusion of the deep IR obser- 
vations obtained with the Spitzer telescope represent a useful 
constraint for the estimate of the physical properties of galaxies 
at high redshift. Last but not least the extensive spectroscopic 
follow up obtained in this field provides a wide set of spectro- 
scopic redshifts. From this public data set we have obtained a 
multi-colour catalogue of galaxies we named GO ODS-MUSIC 
(GOOD S Multicolour South Infrared Catalog, ICrazian et aL 
2006al) . This catalog, where galaxies are selected both in the 



Faber et al. 2005; WiUmer et al 



and other sur- 
20031 



200^ 



Z850 and Ks bands, contains information in 14 bands from the 
U to the Spitzer 8 jim band, and all the available spectroscopic 
information. For all the objects without spectroscopic informa- 
tion we have obtained well calibrated pho tometric re dshifts by 
means of our photometric redshift code (F ontana et al..,200Q) . 

The GOODS-MUSIC catalogue has been already used to 
investigate the physical and clustering properties of high red- 
shift galaxies jpontana e t alJ bOOdlbrazian et al1l2006h. 2007 ; 
Pentericci et al.ll2007l; ICastellano et al.ll2007h . Wifli flie present 
paper we study the galaxy LFs of the red and blue populations, 
enlightening evolutionary features which are characteristic of 
the two populations. 

The paper is organised as follows: in section 2, we describe 
the basic features of our dataset. In section 3, we describe the 
bimodality properties of the sample and we define the loci of 



minimum for the selection of the red/early and the blue/late 
sub-samples as a function of z. In section 4, we compute the 
shape and evolutionary properties of the LFs and the luminosity 
density of both populations. Section 5 is devoted to the analysis 
of the physical properties of the faint early population. 

All the magnitudes used in the present paper are in the AB 
system. An Qa = 0.7, Qm = 0.3, and //q = 70 km s"' Mpc"' 
cosmology is adopted. 

2. THE GOODS-MUSIC SAMPLE 

We use the multicolour catalogue extracted from the southern 
field of the GOODS survey, located in the Chandra Deep Field 
South. The proced ure adopted to extrac t the catalogue is de- 
scribed in detail in^ Grazian et al. ( 2006a ). Here we summarise 
the general features. 

The photometric catalogue was obtained combining 14 im- 
ages from the U-band up to 8 ^m. More specifically it includes 
two U band images from the ESO 2.2 m telescope, an U im- 
age from VLT-VIMOS, the ACS-HST images in four bands 
B, y, / and Z850, the VLT-ISAAC J,H and Ks bands and the 
Spitzer bands at 3.6, 4.5, 5.8 and 8 ^m. All the images anal- 
ysed have an area of 143.2 arcmin^, except for the U-VIMOS 
image (90.2 arcmin^) and the H image (78.0 arcmin^). The 
multicolour catalogue contains 14847 objects, selected either 
in the z and/or in the K s band (version (1.0)). As in previous 



papers (IPoli et al.l 120031 G05) we select galaxies in different 
bands depending on the redshift interval; more specifically we 
select galaxies in the z band at low redshifts (0.2-1.1) and in 
the Ks band at higher redshifts (1.1-3.5). This allows us, as 
explained below, to estimate the galaxy luminosity function in 
the rest fram e 4400 A in the overall red shift inte rval (0.2-3.5) 



As stated bv ICameron & Driven ( 120071) (see also iTruiillo et al 



[2OO6) a careful analysis of the selection efifects due to the detec 
tion c ompleteness is needed. This issue is discussed in the pa- 
per bv lGrazian et alj (l2006al) . In that paper we have evaluated, 
using simulations in the zsso and Ks bands, a 90% complete- 
ness level for elliptical and spiral galaxies of different half-light 
radii and bulge/disk ratios. Since the depth of the image used 
for the galaxy selection varies across the area, a single magni- 
tude limit cannot be defined in each band. As a consequence we 
have divided the z-selected sample and the /T^-selected sam- 
ple in six independent catalogs, each with a well defined area 
and magnitude limit, relative to a 90% completeness level. The 
z-selected catalogs have magnitude limits in the range 24.65- 
26.18, while the magnitude limits in /Ti-selected sample range 
from 21.60 to 23.80, but for most of the sample the typical 
magnitudes limits are zgso ~ 26.18 and Ks ~ 23.5. 

In summary, the z-selected sample has 9862 (after remov- 
ing AGNs and galactic stars) galaxies with about 10% having 
spectroscopic redshift, while the /Ti-selected sample has 293 1 
galaxies with about 27% having spectroscopic redshifts. For 
the galaxies without spectroscopic redshift we use the photo- 
metric on e. Our photometric re dshi ft technique has been de- 
scribed in lGiallongo et alj (Il998l) and lFontana et all ( I2OO0I) . We 
adopt a standard minimisation over a large set of templates 
obtained from synthetic spectral m odels; in particular we use 
those obtained with PEGASE 2.0 dFioc & Rocca-Volmerange 
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1997h as described in iGrazian et alj (|2006aj). The comparison 
with the spectroscopic subsample shows that the accuracy of 
the photometric estimation is very good, with < |Az/(l +z)\ >= 
0.045 in the redshift interval < z < 6. 



As in IPoli et all hoOSi) and lOiallongo et"al] (l2005k great 
care was given to the selection of the sample suited for the esti- 
mate of the Luminosity Function. Indeed we used the z-selected 
sample for galaxies with z < 1 where the 4400 A rest-frame 
wavelength is within or shorter than the zgso-band. For this rea- 
son we included in our LF only galaxies with m[4400(l + z)] < 
ZAB(liin)- This selection guarantees a completeness of the LF 
sample at z < 1 independently of the galaxy colour although 
some galaxies from the original z-limited sample are excluded 
since they have a red spectrum, and consequently a magnitude 
m[4400(l + z)] fainter than our adopted threshold. The same 
procedure was adopted at higher (z=1.0-3.5)redshifts using the 
/Ti-selected sample. The sample selected as described above 
was adopted for all the analysis presented in this paper. 

The method adopted to estimate the rest-frame magnitude 
and the other physical parameters is described in previous 



paper s jPoU et al.l l2003t iGiallongo et al.1 120051: iFontana et al 



2006h . Briefly, it is based on a set of templates, computed wit h 
standard spectral synthesis models ( Bruzual & Charlol 2003 ). 
chosen to broadly encompass the variety of star-formation his- 
tories, metallicities and extinctions of real galaxies. To provide 
a comparison with previous works, we have used the Salpeter 
IMF, ranging over a set of metallicities (from Z = 0.022© 
to Z = 2.5Z0) and dust extinctions (0 < E{B - V) < 1.1, 
wi th a Calzetti extincti on curve). Details are given in Table 1 
of iFontana et al.l (120041) . For each model of this grid, we have 
computed the expected magnitudes in our filter set, and found 
the best-fitting template with a standard minimisation, fix- 
ing the redshift to the spectroscopic or to the photometric one. 
The best-fit parameters of the galaxy were found after scaling 
to the actual luminosity of each observed galaxy. Uncertainties 
in this procedure produced, on average, small errors 10% ) 
in the rest-frame luminosity dElhslllQQTl : iPozzetti et alj|2003b . 
Moreover, the inclusion of the 4 Spitzer bands, longward of 2.2 
//m,for galaxies at z > 2, was essential to sample the spectral 
distribution in the rest-frame optical and near-IR bands, and to 
provide reliable constraints on t he stellar mass and d ust estima- 
tion (for a detailed analysis see Fontana et alj|2006l) . 



3. The semi-analytical model 

In order to make a comparison with current hierarchical mod- 
els of galaxy formation we used our semi-analytical model 
(SAM), described here briefly (for a detailed description see 
Menci et al.lEooil2006h . 

The model connects i) the processes related to the gas 
physics (emission, radiative cooling, heating), ii) the star for- 
mation activity (whose rate is assumed to proceed from the con- 
version of the cold gas mass on a timescale proportional to the 
disk dynamical timescale) and iii) the consequent Supernovae 
heating of the gas to the merging histories of dark matter 
haloes. The model also includes the effect of starbursts trig- 
gered by galaxy interactions and the accretion onto supermas- 
sive black holes at the centre of galaxies powering the AGN 




Fig. 1. left panel: colour-magnitude diagrams in various red- 
shift intervals; the lines represent the best-fit relations for the 
blue and red populations and the locus of the minimum, the 
shaded area is the uncertainty on the minimum. Right panel: 
the histograms of colour distribution projected at Mb - -20 
along the best-fit correlations, the continuous horizontal lines 
are the colour separation at this magnitude, and the dash hori- 
zontal lines are the loci of the maxima. 



activity, with the corresponding energy feedback onto the in- 
terstellar medium. 

We adopt the same choice for the model free parame- 
ters (the normalisation of the star formation timescale and of 
the Supernovae energy feedback) as in the above papers; the 
only changes concern the use of merger trees with a larger 
mass resolution (M - lO'^M©) for progenitors of large mass 
(M > 10'^ Mq) haloes, and the complete depletion of gas in 
haloes with a virial temperature lower t han 4 ■ 10^ K, due to the 
eff'ect of the UV background (see also ISomerville & Primack 
1999h . 



4. BIMODALITY 

4.1. Colour and SSFR properties of the 
GOODS-MUSIC sample 

The recent analysis of the spectral properties of galax- 
ies selected in large or deep surveys has shown the pres- 
ence of a strong bimodality in their colour distribution (e.g. 
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MB(AB)g-19M3(AB)>-lc 




-24 -22 -20 -18 -16 -14 0.4 0.1 
Mb(AB) 

Fig. 2. as in Fig.[T]but for the SSFR-magnitude distribution. 



Strateva et alJ I2OOII: iBaldry et alJ 12004 IWillmer et aP l2006h 



allowing the identification of two main populations, red/early 
and blue/late galaxies mainly on the basis of a single colour, 
e.g. the rest-frame U-V. The local distribution h as bee n stud- 
ied by IStrateva et al] (1200 lb and iBaldrv et al.1 (|2004 in the 



redshifts, by several authors (Bell et al.l 2004; Gialloneo et al. 


2005 


Weiner et al."2005l; Cirasuolo et al. 


2006; .Cucciati et al. 


2006 


Franceschini et al. 


2006). Some effort has been devoted 



in explaining the observed bimodal ity in the framework of 



the hierarchical clusterin g picture ( Menci et al.l 12005 



2006t 



Dekel & BirnboimI I2OO6I) . In particular, Men ci et al.r(l2005h 
proposed that the colour bimodality arises because of two nat- 
ural features: the star formation histories of the massive red 
galaxies, which are formed in biased high-density regions, are 
peaked at higher z as compared to lower mass galaxies; and the 
existence of a non-gravitational mass scale (mo). For m < mo 
the star formation is self regulated and the cold gas content is 
continuously depleted by SN feedback, for m > mo the cold gas 
is not effectively reheated and so the rapid cooling takes place 
at high-z. These different evolu tionary paths led to the present 
day red and blue populations ( Menci et al.ll2005h. When the 
energ y injection from AGN feedback is included ( Menci et al. l 
20061), the bimodal distribution appears at even higher redshifts 



(z > 2). 

Using the rest-frame colour we can separate the red pop- 
ulation from the blue to analyse the evolution of the LFs. A 



recent analysis of the morphological structure of a fraction of 
the GOODS sample has shown a good correlation between the 
red colour and the spheroidal morp hology of galaxies up to 
z ~ 1.5 (see lFranceschini et al ] l2006h . Moreover, as in G05 we 
are interested in the galaxy evolution as a function of the star- 
formation activity. In this respect, the use of the colour criterion 
introduces some population mixing for the red galaxies since it 
is not possible to distinguish an early-type galaxy from a dusty 
star-burst using only the rest- frame t/ - V colour. Therefore, 



we have used the lBruzual & Charloli (l2003h spectral best fit of 
the individual galaxy SEDs to derive the specific star formation 
rate (SSFR) m,/m, (as described in the previous section). We 
are aware that the absolute values in the m,/m, distribution are 
subject to uncertainties due, for example, to the estimate of the 
dust attenuation which depends on the extinction curve adopted 
and to the methods adopted for the mass estimate. We refer to 
our previous paper G05 and references therein for a description 
of the method used and its reliability. However, although some 
degeneracy still remains, we can use this property to separate 
our sample, removing the obvious star-burst galaxies from the 
locus of early-type. 

The results about the colour bimodality at high redshift 
from G05 are here confirmed at a higher statistical level. The 
colour-luminosity relation is shown in Fig.[Tl while the analo- 
gous distribution in SSFR is shown in Fig.|2l 

The minima in the colour-magnitude distribution and in the 
SSFR -magnitude distribution are used to divide the sample in 
red/early and blue/late populations. For an evaluation of this re- 
lation we have adopted the same procedure as in G05. We have 
fit the distribution shown in Fig. [1] and |2] with the sum of two 
gaussians whose mean is a linear function of the absolute mag- 
nitude Mb- Each gaussian has a constant dispersion and each 
sub-sample of galaxies, with a different magnitude limit, has 
been weighted with its covering sky area. Since the statistics 
of the red population is still poor, we have adopted, as in G05, 
the same dependence on the absolute magnitude for both the 
blue and the red populations. Taking into account the different 
normalisations of the two gaussian distributions we have then 
derived the locus of the formal minimum in the sum of the two 
gaussians, separating in this way the two populations. The loci 
of the red/early and blue/late populations are shown in Tab. [1] 
and the minima in Tab.|2] 

The resulting colour distribution projected at Mg - -20 
along the best-fit correlation is also displayed in Fig.[Tl with the 
vertical line indicating the colour separation at that magnitude. 

The same is shown for the SSFR distribution in Fig. |2] 
where the relative numbers of early and late type galaxies can 
also be derived in two different ranges of absolute magnitudes. 

The uncertainty associated with the selection of the min- 
ima has been derived reproducing 100 colour-magnitude plots 
with a MonteCarlo analysis using 100 galaxy catalogs. In each 
catalogue we assigned to each object a different redshift drawn 
from its probability distribution and we associated their rest- 
frame absolute magnitudes and SSFRs. The z probability distri- 
bution naturally takes into account the photometric errors and 
the model degeneracy in the spectral libraries. The uncertainty 
region is shown in figs.[ri|2]as a shaded area, its value is < 0.1 
for the minima in colour, and ~ 0.2 for those in SSFR. 
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The colour/SSFR-magnitude relations for the loci of the 
maxima and minima follow the linear relations < U - V >= 
a ■ (Mb + 20)+ <U-V>2o and <SSFR>^a- (Mb + 20)+ < 
SSFR >2o, whose parameters are listed in Tab.[T]and|2] 

In the colour-magnitude relation we confirm the weak in- 
trinsic blueing with increasing redshift from z ~ 0.4 to z ~ 2 
already found by G05 for both populations although formally 
only at 2cr level. In the redshift bins where the statistics is poor, 
the minimum is extrapolated from the other redshift bins. From 
the bins z = 0.4 - 0.7 and z = 0.7 - 1.0 we extrapolate the 
minimum value < U - V >20- 1.58 in the redshift interval 
Z = 0.2 - 0.4, and from the bins z = 0.7 - 1 .0 and z = 1 .0 - 2.0 
we extrapolate < U - V >2()= 1.23 in the interval z = 2 - 3.5. 

We found no appreciable redshift evolution in the SSFR 
distribution, so in order to increase statistics we have performed 
the fit in the larger redshift intervals z = 0.2 - 1 and 1 - 3.5. 
Moreover, there is no appreciable dependence of SSFR on the 
absolute magnitude, so the colour-magnitude relation is not re- 
lated to similar trends in the specific star-formation rate. 
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4.2. Discussion on tine selection criteria 

One notes in the colour/SSFR-magnitude relations the pres- 
ence of a conspicuous number of intrinsically faint galaxies 
with relatively red colours. They are red with respect to the 
locus of separation of the two populations although, because of 
the colour-magnitude relation, their colours are typical of the 
bright (Mb ~ -22) blue galaxies. In terms of SSFR these galax- 
ies show intermediate values between star forming and early 
type galaxies. The presence of a large number of galaxies be- 
longing to this intermediate population dominates the shape of 
the LF of the red/early type galaxies at the faint end, as shown 
in the next sections. 

However, since the colour dispersion of the blue sequence 
broadens at faint magnitudes, the assumption of a linear pa- 
rameterisation for the minimum could imply a bias for the se- 
lection of the red samp le (for a detailed analysis of the blue 
sequence properties see 



Labbe et al. 2007). For this reason we 



have also performed an analysis deriving the minimum with- 
out any assumption on its dependence on the rest-frame lu- 
minosity. We have concentrated our analysis in the redshift 
interval 0.4-0.7 which have sufficient statistics at faint mag- 
nitudes. We have obtained a volume corrected colour distri- 
bution in four magnitude bins, and for each colour distribu- 
tion we have fit a double-gaussian function as shown in Fig. 
[3] We have adopted the intersection of the two gaussians as 
minimum. In this way we have verified that the locus of the 
minimum is well described by a linear behaviour: performing 
a linear fit to these points we have found the following rela- 
tion < t/ - y >= -0.07 ■ (Mb + 20) + 1.57 (see Fig.©. This 
relation is very similar to our standard analysis and does not 
produce a substantial variation of our results (see also Section 
15.31 and Fig.|8]l. As described above, this analysis can not ex- 
clude a contamination from the star-burst galaxies reddened by 
dust, for this reason we have also performed an analysis of the 
LF on a homogenous sample of galaxies selected using SSFR 
distribution. 



Fig. 3. Volume coiTected colour histogram, in four bins of mag- 
nitude, in the redshift interval 0.4 < z < 0.7. The continuous 
lines are the double-gaussian fit to the colour distributions. 




0.4<redshift<0.7 



-24 



-22 



-20 



-16 



-14 



M,(AB) 



Fig. 4. Points are the minima obtain by the fits shown in Fig [3] 
Continuous line is the linear fit to these minima. 



5. LUMINOSITY FUNCTION 

5.1. The Statistical Analysis 

The LF is computed with the procedure described in G05. We 
have applied to our sample an exten ded version of the standard 
l/Vmax algorithm (.Avni & Bahcalliil980) . As in the previous 
paper, we have used a combination of data derived from regions 
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Table 1. Parameters of the relation between the loci of the max- 
ima and the absolute B magnitude. 



(z) 


a 


^red/early 




U-V 


0.4 0.7 
0.7 1.1 
1.1 2.0 


-0.08±0.01 
-0.08±0.01 
-0.08 


1.79±0.05 
1.74+0.05 
1.69+0.05 


1.09±0.03 
0.93±0.03 
0.90±0.03 


SSFR 


0.2 1.1 
1.1 3.5 


0.019±0.03 
0.019 


-11.59±0.20 
-11.55±0.20 


-9.±0.03 
-9.07±0.03 



Table 2. Parameters of the relation between the locus of the 
minimum and the absolute B magnitude. 





a 


minimum 


U-V 


0.2 0.4 


-0.08 


1.58±0.10* 


0.4 0.7 


-0.08± 0.01 


1.51± 0.07 


0.7 1.1 


-0.08± 0.01 


1.43± 0.06 


1.1 2.0 


-0.08 


1.36± 0.07 


2.0 3.5 


-0.08 


1.23±0.10* 


SSFR 


0.2 1.1 


0.019±0.020 


-10.41+ 0.20 


1.1 3.5 


0.019 


-10.43+ 0.20 



* Extrapolated value 

in the field with different magnitude limits. Indeed, for each 
object and for each j-th region under analysis a set of effec- 
tive volumes Vmaxij) is computed. For a given redshift interval 
(zi,Z2), these volumes are enclosed between zi and Zupij), the 
latter being defined as the minimum between zi and the max- 
imum redshift at which the object could have been observed 
within the magnitude limit of the j-th field. The galaxy number 
density 4>{M, z) in each (Az, AM) bin can then be obtained as: 



/=1 



(1) 



where u){j) is the area in units of steradians corresponding to 
the field j, n is the number of objects in the chosen bin and 
dV/dz is the comoving volume element per steradian. 

The Poisson error in each LF magnitude bin was computed 
adopting the recipe by Gehrelsl (.1986) . valid also for small 
numbers. The uncertainties in the LF due to the photometric 
uncertainties and to the degeneracy of the spectral models used 
to derive redshifts were computed with the same Monte Carlo 
analysis described in the previous section. The uncertainties 
obtained and the Poisson errors have been added in quadrature. 

The l/y,„n.v estimator for the LF can in principle be af- 
fected by small scale galaxy clustering. For this reason a para- 
metric maximum likelihood estimator is also adopted which is 
known to be les s biased respect to small scale clustering (see 
Hevl et al.l[l997l) . 



The parametric analysis of the galaxy LF has been obtained 
from the maximum likelihood analysis assuming for different 
galaxy populations different parameterisation for the LF. The 
maximum likelihood method used here represents an exten- 
sion of the standard ISandage et al.l (Il979h method where sev- 
eral samples can be jointly analysed and where the LF is al- 
lowed to vary with redshift. A more detailed description can be 
found in G05 and a formal derivation of the maximum likeli- 



hood equation is shown in lHevl et al.l (Il997h . 

To describe the B-band LF of the total sample and that of 
blue/late galaxies we assume a Schechter parameterisation: 

cf>(M,z) = 0.4 ■ ln(lO)0*[lO°'^('^*-*')]'+"- 
exp[-10"-4(^'-^)] 

As previous analysis have shown that the evolution of the 
global galaxy LF is manly driven by luminosity evolution, 
we have adopted an M* that is redshi f t dependent (see asl o 



Hevl et al.lll997l: iGiallongo et all 120051: iGabasch et al.ll2006D . 



We have adopted the parameterisation that better matches our 
data in this range of redshifts: 



M*(z) = 



Mo - Ml • z if z < Zcut 
Mo - Ml • Zcut if z > Zcu 



(2) 



The slope a is kept constant with redshift since it is well con- 
strained only at low-intermediate redshift (z < 1). For all the 
sample here considered, we checked that the density parameter 
ip had no significative variation if it had been evolved with the 
parametric form described in G05, thus we keep it as constant. 

As we will see in Fig. |7] and |9] this is not a good descrip- 
tion for the red/early population for which we assume a double 
Schechter form, as fr equently adopted in similar cases in liter- 
ature (e.g. lPopesso et al 2006) : 



cpiM) = ^f(M) + MM) = 
(P* -OA-lnilOy 
([10°-4(*^}-^)]i+«/exp[-10"'*(^/-^'] + 
([100-4(*';-'W)]i+«"exp[-100-4(^J-'^>]) 



(3) 



For a quantitative evaluation of the density evolution at the 
bright end, we have evolved the bright Schechter with a red- 
shift dependent normaUsation: 



<^(M,z) = (pf(M) + norm(z) ■ (t>b{M) 
where 

f 1 if Z < Zcut 

■ iZcut) if Z > Zcut 



(4) 



(5) 



We used the M INUIT package of the CERN library 
("james & Roos 1995 ) for the minimisation. The errors in Tab. 
|3]and|4]are calculated for each parameter, independently of the 
others. 



5.2. B-band luminosity function from z ~ 0.2 to z ~ 3.5 
and the comparison with the other surveys 

The evolution of the total LF is shown in Fig.|5] To compare 
our high z results and our fits with local values we also show 



S. Salimbeni et al.: GOODS-MUSIC: The red and blue galaxy populations 



7 



the local L F from Two-Degree Field Galaxy Redshift Survey 
(2dFGRS; iNorberg et all2()02l dotted line in Fig.|5]l. 

The 1/V,„ax analysis shows that the main kind of evolu- 
tion is due to a brightening of the LF with redshift. We have 
applied the ML analysis to the sample using the evolutionary 
form of the Schechter LF described in eq. |2] where pure lumi- 
nosity evolution is allowed up to a maximum redshift beyond 
which the LF keeps constant. The results in Tab. [3] imply that 
the LF is subject to a mild luminosity evolution only up to z ~ 1 
(AM* ~ 0.7 in the z = 0.2 - 1 interval). At higher z the LF ap- 
pears constant with redshift although at z ~ 3, in the brightest 
bin, a slight excess is present. In any case, the adopted evolu- 
tionary model is acceptable at 2cr level using the standard x~ 
test. 

We also sh ow a comparison with the DEEP2 R edshift 
Survey (D2RS, lOavis etal.ll200l IWillmer et ai]|2006l empty 
points in Fig.|5]l. Although the comparison of the LFs was per- 
formed on data taken from surveys having different magnitude 
limits and redshift estimates (photometric and spectroscopic), 
the agreement is very good in the overlapping magnitude re- 
gions and in all the redshift bins. A general good agreement 
is also found with the LF derived from the VVDS survey by 
Ilbert et al.l (120051) . although from our ML analysis we do not 
have any evidence of steepening of the faint end slope up to 
z ~ 2, as suggested by them. We have compared our results 
with othe r photometric re dshift sur veys like the CO MBO- 17 
survey by Bell et al. (2004) (see alsol Faber et al.l200 5). A good 
agreement is found in the overall redshift interval and in the ap- 
propriate magnitu de interval. The com parison with the FORS 
Deep Field (FDF) dOabasch et al.l2004 ') is less straightforward, 
because of the different redshift intervals used. An acceptable 
agreement is found up to z ~ 1 . At higher redshifts the FDF lu- 
minosity function shows an excess of very bright objects with 
respect to our values as well as the COMBO 17 and DEEP2 
LFs. 

In Fig.|5]we included the LFs derived from our hierarchical 
model described in sec. [3] (see also .Menci et al..,2005, ,2006 ). 
The effect of the UV background is effective in flattening the 
predicted shape of the LF at the faint magn itudes providing an 
agreement better than that obtained before (iPoli et al ] l2003l) . 



5.3. Luminosity function for tlie blue/late and red/early 
galaxies 

In this section we show the shapes and evolutionary be- 
haviours of the luminosity functions derived for the blue/late 
and red/early galaxy populations respectively. We have adopted 
the empirical colour/SSFR selection described in section |4] to 
separate the two populations. 

The shape and redshift evolution of the blue LF is shown in 
Fig.|6]where both the 1 /y„K,i data points and the curves derived 
from the ML analysis are represented. The best fit parameters 
together with their uncertainties are shown in Tab. [3] The LFs 
of the late populations are very similar to the blue ones and we 
do not show the figure. 

As for the total sample, we found that the blue population 
is well represented by the same type of luminosity evolution 




-24 -21 -1 
M, (AB) 



15-24 -21 -If 
Mb (AB) 



Fig. 5. Total LF as a function of redshift. The continuous curves 
come from our maximum likelihood analysis. The dotted line 
is the local LF dNorberg et al. 2002). The filled circles are the 
points obtained with 1 / Vmax method. The em pty circles come 
from the DEEP2 survey (IWillmer et al ]l2006h . The results from 
the COMBO 17 and VDDS surveys are consistent with the 
DEEP2 results and are omitted in the plot. The dashed-point 
line is the model of Menci et al. I (l2006h . 



although faster, with AM* ~ 1.15 in the z - 0.2 - 1 interval. 
The faint end slope appears steeper. This is due of course to the 
fact that the blue population dominates the volume density of 
the total sample at any redshift. 

We first co mpare our results with those of the spec troscopic 
survey D2RS dFaber et al.ll2005l: IWillmer et al.ll2006l) . We note 
that the colour selection they use to separate the two popula- 
tions is based on a f/ - B vs My colour-magnitude relation. We 
have verified in our sample that their colour selection is very 
similar to our criterion. In fact, if we adopt their selection on 
our sample, we reproduce almost the same blue/red galaxy sub- 
samples obtained with our cut. Their LFs are in good agrement 
with our results as shown in Fig. |6] (data with error bars). We 
then compare our results to those of COMBO-17. They use as 
a selection criterion the c olour U -V vs My, which nearly cor- 
responds to that used by IWillmer et al.l d2006t) . The agreement 
with our results is good. 

A direct compar ison with the blue/red LFs by 
Marchesini et al.l d2007h in the redshift interval 2 < z < 3.5 



is not possible since they use a colour separation bluer 
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Table 3. Luminosity Function Parameters for the total and blue/late type sample 



Type 


z 






a 


Ml 


Zcut 




A' 


all 


0.2-3.5 


-20.54+ 


0.12 


-1.34±0.01 


0.93±0.13 


0.96±0.02 


0.0031 


5115 


blue 
late 


0.2-3.5 
0.2-3.5 


-19.97+ 
-20.06+ 


0.20 
0.12 


-1.39±0.02 
-1.36±0.02 


1.43±0.27 
1.29±0.12 


0.98±0.05 
0.99±0.02 


0.0024 
0.0029 


4127 
4612 



than our criterion by 0.2 mag providing a LF with a lower 
normalisation. 

The LF predicted by our hierarchical model was not in- 
cluded in Fig. |6] since it is not appreciably different from that 
of the total sample. 
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Fig. 6. LFs of the blue galaxies as function of redshift. The con- 
tinuous curves comes from our maximum likelihood analysis. 
The dotted line is our fit at z ~ 0.3, reported for comparison 
in all the redshift bins. The big filled circles are the points ob- 
tained with I IV MAX method. The littl e empty circle are points 
from the LFs by lwillmer eTaP (l2006h . 



Concerning the red/early populations our GOODS-MUSIC 

sample allows a sampling of the red LF down to M 16 up 

to z ~ 0.8 providing for the first time a direct evaluation of the 
faint shape of the LF. This is due to the deeper K magnitude 
limit respect to that used in G05, allowing the evaluation of the 
rest frame U-V colour at magnitudes as faint as Mb 16. At 
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Fig. 7. LFs of the red galaxies as function of redshift. The con- 
tinuous curves comes from our maximum likelihood analy- 
sis, the first bin of redshift, which have to low statistic, has 
been excluded from this evolutive analysis. The dotted line 
is our fit at z ~ 0.5, reported for comparison in all the red- 
shift bins. The big filled circles are the points obtained with 
1 I Vmax method. The litt le empty circle are points from the LFs 
by IWilkner et al.1 (l2006h. The dashed-point orange line is the 
model of Menci et al. 



variance with previous works which involve shallower samples, 
a peculiar LF shape is present at z < 0.8 with a minimum and a 
clear upturn at > -18 (Fig.|7]l. 

This overabundance of faint objects with respect to the ex- 
trapolation of the Schechter function was already found in the 
LF of local early type galaxies derived from the 2dF survey 
(Mad gwick et al.l 120021) . Although they used a different sepa- 
ration criterion, their subsample called type 1 is not so differ- 
ent from our subsample, being equivalent to a morphological 
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(AB) 

Fig. 8. LFs of the red galaxies. With the filled circle is indi- 
cated the LF obtained with the linear color selection in Tab. 
|2] The empty circles indicate the LF obtained from the color 
histograms in bins of magnitude (see Fig.[3]and Fig.HJi. 



sample of E, SO and Sa. A similar upturn was also found in 
the local LF of red g alaxies derived from the Sloan survey by 
Blanton et all ilOO^ . 

We have checked that the characteristic shape found is 
not critically dependent on the specific choice of the colour- 
magnitude or SSFR -magnitude separation. Indeed, changing 
the parameters of the linear fit in colour-SSFR separation at 
2-cr level, the shape and, in particular, the upturn do not change 
appreciably. We have also found that the overall shape of the 
LF does not change, if we use a separation between blue and 
red galaxies found with an analysis without any assumption 
about the parametric dependence on the the rest-frame magni- 
tude (see section |4~2] |. Moreover, fitting with a double gaussian 
function the bimodal distribution in the faintest magnitude in- 
terval (see Fig. [3] -17 < MB < -15) the contamination by 
any blue population in the locus of the red population is not 
so strong. In other words the FWHM of the blue gaussian dis- 
tribution is relatively narrow. In this case the expected fraction 
of blue galaxies expected redward of the selected minimum is 
only 14% of the red population in the same colour region. For 
this reason the shape of the red LF shown in Fig. 8 relative to 
the redshift interval 0.4 < z < 0.6 remains almost unchanged in 
the two fainter magnitude bins if we remove this small fraction. 

Also for the red population we have compared our LFs with 
that derived from the major surveys of colour selected galax- 
ie s. In Fig. |7l we show the LF from the spectroscopic survey 
of Willmer et al. For the red galaxies the agreement is 

good for Mb < -20, wh ere the incom pleteness is negligible in 
their sample (see fig 8, IWillmer et al.i,2006) . Their shallower 
sample can not probe the raise at the faint end present in our 



-5 
-2 

-3 

-4 

-5 
-2 



1 1 1 1 1 1 1 1 1 1 
; 0.2<z<0.4 ^ 

; : 
" 1 , , 1 , , 1 , , 1 ■ 


1 1 1 1 1 1 1 1 1 1 
; 0.4<z<0.6 ^JJH ^ 

: 1 '1 1 1 1 1 1 1 1 1 " 


_ 

: 0- 


1 1 1 1 1 1 1 
6<z<0.8 ■ 


■ 0.8<z<l.l 


1 i 1 1 i 1 1 1 1 1 

\ 1.1<Z<2 

■ 1 ■'/ 


: Ml 1 1 1 1 1 1 1 1 " 

'■ 2<z<3.5 

" 1 ■■/ , 1 , , 1 , , 1 " 



-24 -21 -18 -15 -24 -21 -18 
M„ (AB) M„ (AB) 



-15 



Fig. 9. LFs of the early type galaxies as function of redshift. 
The continuous curves comes from our maximum likelihood 
analysis. The dotted line is our fit at z ~ 0.5, reported for com- 
parison in all the redshift bins. The big filled circles are the 
points obtained with 1 /Vmax method. 

de eper data. The sa me holds for the two ph otometric surveys 
of lBell et alj (|2004 and lBrown etakl (|2007|) . 

Concerning the parametric analysis of the evolutionary LFs 
of the red/early galaxy population, given the excess of faint ob- 
jects a single Schechter shape does not provide an acceptable 
description of the data. For this reason we have adopted a dou- 
ble Schechter function as described in eq. H] The best fit pa- 
rameters are shown in Tab. |4] The best fit value of the brighter 
Schechter slope is rather flat (ai, ~ -0.7) in agreement with 
what found in G05 and in the Bell et al. sample. The fainter 
slope is steeper approaching the value a/ ~ -1.8. As for the 
redshift evolution we have adopted the density evolution law 
described in eq.|5] where the Schechter shape at the faint end is 
kept constant at all redshifts. The brighter one is constant only 
up to a given redshift z„„ beyond which decreases as a power 
law in redshift. We find a constant density up to z ~ 0.7 and 
thereafter a decrease by a factor ~ 5 up to z ~ 3.5. 

The LF evolution of the early galaxies selected from their 
SSFR value is very similar to that of the red ones although the 
high redshift density evolution is more pronounced with a de- 
crease by a factor ~ 10 in the interval z - 0.7 - 3.5. This differ- 
ence is caused by the presence, in the red sample at higher red- 
shift, of a high fraction of galaxies having SEDs consistent with 
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Table 4. Luminosity Function Parameters for red and early type 
galaxies, in the redshift interval 0.4-3.5 



parameter 


Red 


Early 


N 


913 


472 


<t>l 


0.0020 


0.0017 


7 


-1.045±0.13 


-1.45+0.16 




0.65+0.01 


0.67±0.01 




-0.76±0.06 


-0.40±0.1 


m; 


-21.29±0.1 


-21.22±0.12 


«/ 


-1.77±0.2 


-1.54+0.22 




-17.04±0.12 


-17.06+0.12 



those of a dusty and starburst galaxy. This fraction amounts to 
~ 70% at Mb ~ -21.5 andz ~ 3. 

Finally, the comparison with our hierarchical CDM model 
shows a slight flattening of the LF at intermediate luminosi- 
ties. This is because the red population is mainly contributed 
by galaxies with larger M/L ratios; these are mainly attained 
in massive objects (due to the ineffectiveness of gas cooling, 
to their earlier conversion of gas into stars, and to the effect of 
AGN feedback) or in low-mass objects (due to the gas deple- 
tion originated from the different feedback mechanisms, partic- 
ulary effective in shallow potential wells). However, the model 
still overpredicts the LF at faint magnitudes; we shall investi- 
gate the origin of such an effect (probably originated at high- 
redshifts) in a future paper. 

5.4. Luminosity densities 

To compare in a global way the redshift evolution of the blue 
and red galaxies we have computed the B band luminosity den- 
sities of the two populations as a function of redshift. To make 
the comparison homogeneous we have computed the contribu- 
tion of the same bright population with Mb{AB) < -20.2 at all 
redshifts. The results are shown in Fig.[TO]for the blue/red and 
late/early populations. The redshift bins are selected to have 
a comparable number of objects in the considered magnitude 
range. In fact for the highest redshift bin the lowest luminosity 
of the data is Mb(AB) = -21. For this reason we have added 
the contribution of the fainter sources using the extrapolation 
of the parametric LF. 

The uncertainty associated with the luminosity density is 
the sum in quadrature of the error estimated through jack- 
knifj^Efron .1982.) analysis and of the error obtained from 
the MonteCarlo analysis as described in sec. 15.11 The first 
contribution is associated with the clustering properties of 
the field and the second with the photometric uncertainties. 
Local lu minosity densities, obt ained from the integration of the 
LFs by Madg wick et"ai] (l2002h (z=0.04) and lBell et all (l2003b 
(z=0.07), are also included for comparison. 

The blue/late galaxies increase steadily their luminosity 
density up to z ~ 3.5 while the luminosity density of the 
red/early population is nearly constant up to z ~ 1 and then de- 
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Fig. 10. Upper panel: Galaxy luminosity density for red (filled 
square) and blue galaxies (empty square). Lower panel: the 
same of the upper panel but for early type galaxies (filled 
square) and late type galaxies (empty squar e). Small points 



are the local lumino sity density esti mates by iMadgwick et al 



(l2002h (z=0.04) and lBeUetalJ (120031) (z=0.07) 



creases by a factor ~ 3 at z ~ 3.5. This confirms our previous 
result presented in G05, which shows an appreciable decline of 
the bright red/early population only at z > 1 . As already noted 
this is not in qualitative contrast with the hierarchical scenario, 
since the bright and hence massive early population is devel- 
oping early in the cosmic time in specific overdense regions 
where the evolution is also accelerated by merging processes. 
However, the detailed quantitative agreement of the LFs pre- 
dicted by specific models depends on the details of the main 
physical processes and a satisfactory agreement is not yet ob- 
tained especially at the fainter magnitudes (see previous sec- 
tion). 

6. Red/early faint galaxies 

As shown in the previous section, the presence of an upturn 
in the LF of the red galaxy population at faint magnitudes is a 
new feature emerging from the analysis of deep NIR selected 
galaxies, with respect to previous works at this redshift. To de- 
rive hints on the nature of the population responsible for this 
excess we have analysed their colour and spatial properties. 

The peculiar shape of the LF represented by a double 
Schechter form appears similar to that obtained for galaxies 
in local rich clusters. Indeed recent studies of the total galaxy 
luminosity functions in clusters selected from the RASS-SDSS 
survey show an upturn at fain t Mr which depends on the dis- 
tance from the cluster center ( Popesso et al. ' 2006). The main 



* The jackknife analysis is performed recomputing the statistic es- 
timate leaving out one observation at a time from the sample. 



contribution to this local excess comes from th e red population 
selected with u-r > 2.22 (IStrateva et al .l200lh . They also noted 
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galaxy in the B band among the objects we are considering. At 
z > Zc the fraction of detected objects is: 
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Fig. 11. panel a: Stellar mass distribution for the early and late 
population, the area of each histogram is normalised to unity. 
The arrows indicate the mean value of the stellar mass for the 
early population (thin arrow) and for the late population (tick 
arrow) panel b: as the panel a but for the distribution of pipmed 
panel c: Early type galaxies vs. late type galaxies fraction as 
function of the density contrast (p/pmed)- 



that the ratio between red and blue galaxies increases with the 
density in the clusters. 

It is interesting to explore whether a dependence on the 
environment is present for our faint early subsample. To this 
end we have adopted a 2-D density analysis of the 20 first- 
neighbour method. We describe here briefly t he procedure 



adopte d to derive densities, while we refer to iTrevese et al 
(l2007h "or a detailed description of the method. A 2-D density 
has been assigned to each object and a density map has been 
derived in the field in the redshift interval z = 0.4 - 0.6. 
A surface density 



S„ = n/(nDl„) 



(6) 



of galaxies was computed considering the projected distances 
Dnn to the n-th nearest neighbo ur, as done by Dressier 



( lDressleJl980t iDressler et al .Il997h . We divide the survey area 
in cells whose extension depends on the observational accu- 
racy. For each cell we count neighbouring objects at increasing 
radial distance until a number n = 20 of objects is reached. 
In counting galaxies we must take into account the increase 
of limiting luminosity with increasing redshift for a given flux 
limit. If m/„„ is the limiting (apparent) magnitude in a fixed ob- 
serving band, at each redshift z we detect only objects brighter 
than an absolute magnitude Mumiz), decreasing (brightening) 
with z ■ We assume the lower limit of the interval as reference 
redshift Zc below which we detect all objects brighter than the 
relevant = MumiZc), which is the magnitude of the faintest 



*(z) = 



(7) 



where <I)(M) is the type and z-dependent galaxy B-band lumi- 
nosity function presented above. Thus, in evaluating the galax- 
ies number density, we apply a limiting magnitude correction 
by assigning a weight w{z) = 1 / s(z) to each detected galaxy 
of redshift z- In this way we have evaluated the surface density 
field and assigned a density value to each object. 

For our analysis we have selected galaxies in three magni- 
tude regions; an intermediate region, -19 < Mb{AB) < -17, 
where the LF of the early population is flat, and two external 
steeper regions at the bright and faint end of the LF, MsiAB) < 
-19 and Mb{AB) > -17 respectively. 

First of all we note that early galaxies represent the most 
massive galaxies in each luminosity interval (Fig. [TT]panel a). 
In particular even at the faint end of the LF the early popula- 
tion is clearly segregated in stellar mass with values one order 
of magnitude greater on average with respect to the late popu- 
lation. 

Looking to the brightest fraction, a clear difference as a 
function of the density of the environment is found between the 
early and the late populations. In particular early galaxies tend 
to populate regions of higher density. This is shown in panels 
b,c where the two distributions are represented as a function of 
the density of the environment. The ratio early /late increases 
with density since the average density of the early galaxies is 
somewhat greater (1.4) with respect to the one of the late pop- 
ulation. 

This different behaviour becomes less evident with decreas- 
ing luminosity and almost disappears at the faint end of the two 
LFs. We note in this respect that the limited area covered by our 
sample does not allow an evaluation of the environment depen- 
dence up to the high densities typical of clusters, like those 
probed by e.g. the Sloan survey. 

Thus, the scenario that emerges is one where major evolu- 
tionary differences between the early and late populations act 
in the relatively bright galaxies with M < -17 producing the 
largest differences in the shapes of the two LFs in the inter- 
val -21 < Mb{AB) < -18 (flat shape for the early, steep for 
the late). At faint magnitudes the two LFs tend to converge to 
the same volume density. From this analysis the characteristic 
shape of the red/early LF does not seem to depend strongly on 
the environmental properties. 

7. Summary 

We have used a composite sample of galaxies selected in deep 
NIR images, obtained from the GOODS public survey, to study 
the evolution of the galaxy LFs of red/early and blue/late 
galaxy populations selected using the colour and SSFR statis- 
tical properties of the sample. 

The observed U - V colour and SSFR distributions show 
a clear bimodality up to z ~ 3, confirming the results obtained 
in G05 at a higher level of statistical confidence. We found a 
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trend with redshift for the colour magnitude distribution with 
an intrinsic blueing of about 0.15 mag. in the redshift interval 
Z — 0.4-2.0 for both populations. This observed bimodality can 
be explained in a hierarchical clustering scenario as due to the 
different star for mation histories of the red/early and blue/late 
galaxies, see e.g. Menci et al. ( 20061) . 

For the total and the blue/late sample the LF is well de- 
scribed by a Schechter function and shows a mild luminosity 
evolution in the redshift interval z - 0.2 - 1 (e.g. AM* ~ 0.7 
for the total sample; AM* ~ 1 for the blue/late fraction), 
while at higher redshifts the LFs are consistent with no evolu- 
tion. A comparison with our hierarchical CDM model shows a 
good agreement at bright and intermediate magnitudes. A bet- 
ter agreement of the model has also been found at fainter mag- 
nitudes due to the suppression of star formation in small objects 
by the action of an ionising UV background. 

The shape of the red/early luminosity function is better con- 
strained only at low and intermediate redshifts and it shows an 
excess of faint red dwarfs with respect to the extrapolation of 
a flat Schechter function. In fact a minimum around magnitude 
Mb{AB) - -18 is present together with an upturn at fainter 
magnitudes. This peculiar shape has been represented by the 
sum of two Schechter functions. 

We found that the bright one is constant up to z ~ 0.7 be- 
yond which it decreases in density by a factor ~ 5 (10 for the 
early galaxies) up to redshift z ~ 3.5. The comparison with 
our hierarchical CDM model shows that, although the predicted 
LF has a slight flattening at intermediate luminosity, the model 
still overpredicts the LF at faint magnitudes. The bright end 
of blue and red LFs at low and intermediate redshifts is in 
good agreement with recent estimates from the DEEP2 spec- 
troscopic survey. As a consequence of this complex evolution- 
ary behaviour, the luminosity densities of the relatively bright 
(Mb{AB) < -20.2) red/early and blue/late galaxies show a bi- 
furcation beyond redshift z ~ I - Indeed the LD of the blue/late 
population keeps increasing up to z ~ 3.5, while the luminos- 
ity density of red/early galaxies decreases by a factor ~ 2 - 3 
respectively in the z = 1 - 3.5 interval. 

To derive hints on the nature of the galaxies responsible 
for the peculiar shape of the red/early LF, we have performed 
an analysis of their stellar masses and spatial distribution. We 
found that the early galaxies have systematically higher stel- 
lar masses with respect to the late ones for a given B band 
luminosity. Brighter early galaxies have a spatial distribution 
more concentrated in higher density regions if compared to the 
late ones of the same luminosity class. On the contrary, fainter 
early and late galaxies show a very similar spatial distribution. 
Thus, the different environmental properties do not seem to be 
the main responsible for the difference in shape at intermediate 
magnitudes between the blue and red LFs. The latter seem to 
stem from the different star formation and feedback histories 
corresponding to different possible merging trees (evolutionary 
paths) leading to the final assembled galaxy; this specific his- 
tory, driving the evolution of the star formation, leads to the 
different M/L ratios characterising the different properties of 
blue/late and red/early galaxies. In summary, the peculiar shape 
of the red LF is mainly driven by the nature of the galaxy merg- 
ing tree rather than by the nurture where the galaxy has grown. 
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